  function G= var_comp_one(rho);
% function G= var_comp_one(rho);
% Compute companion form a univariate AR(p) process with RHO lags 
% ================================================================
rho = rho(:); 
np = length( rho ); 
if np == 1 
    G=rho; 
end 
G = zeros( np ); 
G(1,:) = rho'; 
G(2:end,1:end-1) = eye( np-1 ); 
